# P` L ` O ` T ` T ` I ` N ` G

# we are using this for plot: ###### XXXX_-XX_010





#######		#######		#######		#######		#######
#######		#######		#######		#######		#######
#######		#######		#######		#######		#######
#######		#######		#######		#######		#######
#######		#######		#######		R M S E
#######		#######		#######		R M S E
#######		#######		#######		R M S E
#######		#######		#######		R M S E
#######		#######		#######		#######		#######
#######		#######		#######		#######		#######
#######		#######		#######		#######		#######
#######		#######		#######		#######		#######
#######		#######		#######		#######		#######

rm(list=ls())
setwd("/Users/lleemann/Desktop/Replication Files/Simulation/Resulting data")
load("Data_Sample_1000.RData")

objects()

par(family="CMU Serif", mfrow=c(1,3))
#par(family="CMU Serif")
plot(mean(correlation.1000[,6]), sqrt(mean((as.numeric(Param.struc.1000[,6,2])-1)^2)), xlim=c(-1,1), ylim=c(-0.05,0.35), pch=16, col=rgb(255,165,0,250,maxColorValue=255), cex=2.6, ylab="Root Mean Squared Error", xlab="Correlation of Errors", main=expression(paste("Estimate of ", beta[21])))
points(mean(correlation.1000[,6]),sqrt(mean((as.numeric(param.selec.1000[,6,2])-1)^2)),pch=16, col=rgb(0,139,139,200,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,6]),sqrt(mean((as.numeric(param.selStruc.1000[,6,2])-1)^2)),pch=16, col=rgb(160,32,240,200,maxColorValue=255), cex=2.6)
abline(h=0, lty=2)

points(mean(correlation.1000[,7]), sqrt(mean((as.numeric(Param.struc.1000[,6,6])-1)^2)), pch=16, col=rgb(255,165,0,250,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,7]),sqrt(mean((as.numeric(param.selec.1000[,6,6])-1)^2)),pch=16, col=rgb(0,139,139,200,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,7]),sqrt(mean((as.numeric(param.selStruc.1000[,6,6])-1)^2)),pch=16, col=rgb(160,32,240,200,maxColorValue=255), cex=2.6)


points(mean(correlation.1000[,8]), sqrt(mean((as.numeric(Param.struc.1000[,6,10])-1)^2)), pch=16, col=rgb(255,165,0,250,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,8]),sqrt(mean((as.numeric(param.selec.1000[,6,10])-1)^2)),pch=16, col=rgb(0,139,139,200,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,8]),sqrt(mean((as.numeric(param.selStruc.1000[,6,10])-1)^2)),pch=16, col=rgb(160,32,240,200,maxColorValue=255), cex=2.6)


points(mean(correlation.1000[,9]), sqrt(mean((as.numeric(Param.struc.1000[,6,14])-1)^2)), pch=16, col=rgb(255,165,0,250,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,9]),sqrt(mean((as.numeric(param.selec.1000[,6,14])-1)^2)),pch=16, col=rgb(0,139,139,200,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,9]),sqrt(mean((as.numeric(param.selStruc.1000[,6,14])-1)^2)),pch=16, col=rgb(160,32,240,200,maxColorValue=255), cex=2.6)


points(mean(correlation.1000[,10]), sqrt(mean((as.numeric(Param.struc.1000[,6,18])-1)^2)), pch=16, col=rgb(255,165,0,250,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,10]),sqrt(mean((as.numeric(param.selec.1000[,6,18])-1)^2)),pch=16, col=rgb(0,139,139,200,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,10]),sqrt(mean((as.numeric(param.selStruc.1000[,6,18])-1)^2)),pch=16, col=rgb(160,32,240,200,maxColorValue=255), cex=2.6)

#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
# Lines for the plot
# yellow
segments(x0=mean(correlation.1000[,6]), y0=sqrt(mean((as.numeric(Param.struc.1000[,6,2])-1)^2)), x=mean(correlation.1000[,7]), y=sqrt(mean((as.numeric(Param.struc.1000[,6,6])-1)^2)), lty=2, col=rgb(255,165,0,250,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,7]), y0=sqrt(mean((as.numeric(Param.struc.1000[,6,6])-1)^2)), x=mean(correlation.1000[,8]), y=sqrt(mean((as.numeric(Param.struc.1000[,6,10])-1)^2)), lty=2, col=rgb(255,165,0,250,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,8]), y0=sqrt(mean((as.numeric(Param.struc.1000[,6,10])-1)^2)), x=mean(correlation.1000[,9]), y=sqrt(mean((as.numeric(Param.struc.1000[,6,14])-1)^2)), lty=2, col=rgb(255,165,0,250,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,9]), y0=sqrt(mean((as.numeric(Param.struc.1000[,6,14])-1)^2)), x=mean(correlation.1000[,10]), y=sqrt(mean((as.numeric(Param.struc.1000[,6,18])-1)^2)), lty=2, col=rgb(255,165,0,250,maxColorValue=255), lwd=2)

# green
segments(x0=mean(correlation.1000[,6]), y0=sqrt(mean((as.numeric(param.selec.1000[,6,2])-1)^2)), x=mean(correlation.1000[,7]), y=sqrt(mean((as.numeric(param.selec.1000[,6,6])-1)^2)), lty=2, col=rgb(0,139,139,200,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,7]), y0=sqrt(mean((as.numeric(param.selec.1000[,6,6])-1)^2)), x=mean(correlation.1000[,8]), y=sqrt(mean((as.numeric(param.selec.1000[,6,10])-1)^2)), lty=2, col=rgb(0,139,139,200,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,8]), y0=sqrt(mean((as.numeric(param.selec.1000[,6,10])-1)^2)), x=mean(correlation.1000[,9]), y=sqrt(mean((as.numeric(param.selec.1000[,6,14])-1)^2)), lty=2, col=rgb(0,139,139,200,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,9]), y0=sqrt(mean((as.numeric(param.selec.1000[,6,14])-1)^2)), x=mean(correlation.1000[,10]), y=sqrt(mean((as.numeric(param.selec.1000[,6,18])-1)^2)), lty=2, col=rgb(0,139,139,200,maxColorValue=255), lwd=2)

# purple
segments(x0=mean(correlation.1000[,6]), y0=sqrt(mean((as.numeric(param.selStruc.1000[,6,2])-1)^2)), x=mean(correlation.1000[,7]), y=sqrt(mean((as.numeric(param.selStruc.1000[,6,6])-1)^2)), lty=2, col=rgb(160,32,240,200,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,7]), y0=sqrt(mean((as.numeric(param.selStruc.1000[,6,6])-1)^2)), x=mean(correlation.1000[,8]), y=sqrt(mean((as.numeric(param.selStruc.1000[,6,10])-1)^2)), lty=2, col=rgb(160,32,240,200,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,8]), y0=sqrt(mean((as.numeric(param.selStruc.1000[,6,10])-1)^2)), x=mean(correlation.1000[,9]), y=sqrt(mean((as.numeric(param.selStruc.1000[,6,14])-1)^2)), lty=2, col=rgb(160,32,240,200,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,9]), y0=sqrt(mean((as.numeric(param.selStruc.1000[,6,14])-1)^2)), x=mean(correlation.1000[,10]), y=sqrt(mean((as.numeric(param.selStruc.1000[,6,18])-1)^2)), lty=2, col=rgb(160,32,240,200,maxColorValue=255), lwd=2)


legend(.4,.3, legend=c("Stc", "HSl", "StSl"), col=c(rgb(255,165,0,250,maxColorValue=255),rgb(0,139,139,200,maxColorValue=255),rgb(160,32,240,200,maxColorValue=255)), cex=1.2, pch=16, bty="n")


### 2nd plot

plot(mean(correlation.1000[,6]), sqrt(mean((as.numeric(Param.struc.1000[,5,2])-1)^2)), xlim=c(-1,1), ylim=c(-0.05,0.35), pch=16, col=rgb(255,165,0,250,maxColorValue=255), cex=2.6, ylab="Root Mean Squared Error", xlab="Correlation of Errors", main=expression(paste("Estimate of ", beta[22])))
points(mean(correlation.1000[,6]),sqrt(mean((as.numeric(param.selec.1000[,5,2])-1)^2)),pch=16, col=rgb(0,139,139,200,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,6]),sqrt(mean((as.numeric(param.selStruc.1000[,5,2])-1)^2)),pch=16, col=rgb(160,32,240,200,maxColorValue=255), cex=2.6)
abline(h=0, lty=2)

points(mean(correlation.1000[,7]), sqrt(mean((as.numeric(Param.struc.1000[,5,6])-1)^2)), pch=16, col=rgb(255,165,0,250,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,7]),sqrt(mean((as.numeric(param.selec.1000[,5,6])-1)^2)),pch=16, col=rgb(0,139,139,200,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,7]),sqrt(mean((as.numeric(param.selStruc.1000[,5,6])-1)^2)),pch=16, col=rgb(160,32,240,200,maxColorValue=255), cex=2.6)


points(mean(correlation.1000[,8]), sqrt(mean((as.numeric(Param.struc.1000[,5,10])-1)^2)), pch=16, col=rgb(255,165,0,250,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,8]),sqrt(mean((as.numeric(param.selec.1000[,5,10])-1)^2)),pch=16, col=rgb(0,139,139,200,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,8]),sqrt(mean((as.numeric(param.selStruc.1000[,5,10])-1)^2)),pch=16, col=rgb(160,32,240,200,maxColorValue=255), cex=2.6)


points(mean(correlation.1000[,9]), sqrt(mean((as.numeric(Param.struc.1000[,5,14])-1)^2)), pch=16, col=rgb(255,165,0,250,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,9]),sqrt(mean((as.numeric(param.selec.1000[,5,14])-1)^2)),pch=16, col=rgb(0,139,139,200,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,9]),sqrt(mean((as.numeric(param.selStruc.1000[,5,14])-1)^2)),pch=16, col=rgb(160,32,240,200,maxColorValue=255), cex=2.6)


points(mean(correlation.1000[,10]), sqrt(mean((as.numeric(Param.struc.1000[,5,18])-1)^2)), pch=16, col=rgb(255,165,0,250,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,10]),sqrt(mean((as.numeric(param.selec.1000[,5,18])-1)^2)),pch=16, col=rgb(0,139,139,200,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,10]),sqrt(mean((as.numeric(param.selStruc.1000[,5,18])-1)^2)),pch=16, col=rgb(160,32,240,200,maxColorValue=255), cex=2.6)

#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
# Lines for the plot
# yellow
segments(x0=mean(correlation.1000[,6]), y0=sqrt(mean((as.numeric(Param.struc.1000[,5,2])-1)^2)), x=mean(correlation.1000[,7]), y=sqrt(mean((as.numeric(Param.struc.1000[,5,6])-1)^2)), lty=2, col=rgb(255,165,0,250,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,7]), y0=sqrt(mean((as.numeric(Param.struc.1000[,5,6])-1)^2)), x=mean(correlation.1000[,8]), y=sqrt(mean((as.numeric(Param.struc.1000[,5,10])-1)^2)), lty=2, col=rgb(255,165,0,250,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,8]), y0=sqrt(mean((as.numeric(Param.struc.1000[,5,10])-1)^2)), x=mean(correlation.1000[,9]), y=sqrt(mean((as.numeric(Param.struc.1000[,5,14])-1)^2)), lty=2, col=rgb(255,165,0,250,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,9]), y0=sqrt(mean((as.numeric(Param.struc.1000[,5,14])-1)^2)), x=mean(correlation.1000[,10]), y=sqrt(mean((as.numeric(Param.struc.1000[,5,18])-1)^2)), lty=2, col=rgb(255,165,0,250,maxColorValue=255), lwd=2)

# green
segments(x0=mean(correlation.1000[,6]), y0=sqrt(mean((as.numeric(param.selec.1000[,5,2])-1)^2)), x=mean(correlation.1000[,7]), y=sqrt(mean((as.numeric(param.selec.1000[,5,6])-1)^2)), lty=2, col=rgb(0,139,139,200,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,7]), y0=sqrt(mean((as.numeric(param.selec.1000[,5,6])-1)^2)), x=mean(correlation.1000[,8]), y=sqrt(mean((as.numeric(param.selec.1000[,5,10])-1)^2)), lty=2, col=rgb(0,139,139,200,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,8]), y0=sqrt(mean((as.numeric(param.selec.1000[,5,10])-1)^2)), x=mean(correlation.1000[,9]), y=sqrt(mean((as.numeric(param.selec.1000[,5,14])-1)^2)), lty=2, col=rgb(0,139,139,200,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,9]), y0=sqrt(mean((as.numeric(param.selec.1000[,5,14])-1)^2)), x=mean(correlation.1000[,10]), y=sqrt(mean((as.numeric(param.selec.1000[,5,18])-1)^2)), lty=2, col=rgb(0,139,139,200,maxColorValue=255), lwd=2)

# purple
segments(x0=mean(correlation.1000[,6]), y0=sqrt(mean((as.numeric(param.selStruc.1000[,5,2])-1)^2)), x=mean(correlation.1000[,7]), y=sqrt(mean((as.numeric(param.selStruc.1000[,5,6])-1)^2)), lty=2, col=rgb(160,32,240,200,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,7]), y0=sqrt(mean((as.numeric(param.selStruc.1000[,5,6])-1)^2)), x=mean(correlation.1000[,8]), y=sqrt(mean((as.numeric(param.selStruc.1000[,5,10])-1)^2)), lty=2, col=rgb(160,32,240,200,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,8]), y0=sqrt(mean((as.numeric(param.selStruc.1000[,5,10])-1)^2)), x=mean(correlation.1000[,9]), y=sqrt(mean((as.numeric(param.selStruc.1000[,5,14])-1)^2)), lty=2, col=rgb(160,32,240,200,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,9]), y0=sqrt(mean((as.numeric(param.selStruc.1000[,5,14])-1)^2)), x=mean(correlation.1000[,10]), y=sqrt(mean((as.numeric(param.selStruc.1000[,5,18])-1)^2)), lty=2, col=rgb(160,32,240,200,maxColorValue=255), lwd=2)




legend(.4,.3, legend=c("Stc", "HSl", "StSl"), col=c(rgb(255,165,0,250,maxColorValue=255),rgb(0,139,139,200,maxColorValue=255),rgb(160,32,240,200,maxColorValue=255)), cex=1.2, pch=16, bty="n")
####

#### 


####

#### 


####

#par(family="CMU Serif", mfrow=c(1,2))
par(family="CMU Serif")
plot(mean(correlation.1000[,6]), sqrt(mean((as.numeric(Param.struc.1000[,3,2])-1)^2)), xlim=c(-1,1), ylim=c(-0.1,0.7), pch=16, col=rgb(255,165,0,250,maxColorValue=255), cex=2.6, ylab="Root Mean Squared Error", xlab="Correlation of Errors", main=expression(paste("Estimate of ", beta[11])))
points(mean(correlation.1000[,6]),sqrt(mean((as.numeric(param.selec.1000[,3,2])-1)^2)),pch=16, col=rgb(0,139,139,200,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,6]),sqrt(mean((as.numeric(param.selStruc.1000[,3,2])-1)^2)),pch=16, col=rgb(160,32,240,200,maxColorValue=255), cex=2.6)
abline(h=0, lty=2)

points(mean(correlation.1000[,7]), sqrt(mean((as.numeric(Param.struc.1000[,3,6])-1)^2)), pch=16, col=rgb(255,165,0,250,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,7]),sqrt(mean((as.numeric(param.selec.1000[,3,6])-1)^2)),pch=16, col=rgb(0,139,139,200,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,7]),sqrt(mean((as.numeric(param.selStruc.1000[,3,6])-1)^2)),pch=16, col=rgb(160,32,240,200,maxColorValue=255), cex=2.6)


points(mean(correlation.1000[,8]), sqrt(mean((as.numeric(Param.struc.1000[,3,10])-1)^2)), pch=16, col=rgb(255,165,0,250,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,8]),sqrt(mean((as.numeric(param.selec.1000[,3,10])-1)^2)),pch=16, col=rgb(0,139,139,200,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,8]),sqrt(mean((as.numeric(param.selStruc.1000[,3,10])-1)^2)),pch=16, col=rgb(160,32,240,200,maxColorValue=255), cex=2.6)


points(mean(correlation.1000[,9]), sqrt(mean((as.numeric(Param.struc.1000[,3,14])-1)^2)), pch=16, col=rgb(255,165,0,250,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,9]),sqrt(mean((as.numeric(param.selec.1000[,3,14])-1)^2)),pch=16, col=rgb(0,139,139,200,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,9]),sqrt(mean((as.numeric(param.selStruc.1000[,3,14])-1)^2)),pch=16, col=rgb(160,32,240,200,maxColorValue=255), cex=2.6)


points(mean(correlation.1000[,10]), sqrt(mean((as.numeric(Param.struc.1000[,3,18])-1)^2)), pch=16, col=rgb(255,165,0,250,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,10]),sqrt(mean((as.numeric(param.selec.1000[,3,18])-1)^2)),pch=16, col=rgb(0,139,139,200,maxColorValue=255), cex=2.6)
points(mean(correlation.1000[,10]),sqrt(mean((as.numeric(param.selStruc.1000[,3,18])-1)^2)),pch=16, col=rgb(160,32,240,200,maxColorValue=255), cex=2.6)

#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
# Lines for the plot
# yellow
segments(x0=mean(correlation.1000[,6]), y0=sqrt(mean((as.numeric(Param.struc.1000[,3,2])-1)^2)), x=mean(correlation.1000[,7]), y=sqrt(mean((as.numeric(Param.struc.1000[,3,6])-1)^2)), lty=2, col=rgb(255,165,0,250,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,7]), y0=sqrt(mean((as.numeric(Param.struc.1000[,3,6])-1)^2)), x=mean(correlation.1000[,8]), y=sqrt(mean((as.numeric(Param.struc.1000[,3,10])-1)^2)), lty=2, col=rgb(255,165,0,250,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,8]), y0=sqrt(mean((as.numeric(Param.struc.1000[,3,10])-1)^2)), x=mean(correlation.1000[,9]), y=sqrt(mean((as.numeric(Param.struc.1000[,3,14])-1)^2)), lty=2, col=rgb(255,165,0,250,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,9]), y0=sqrt(mean((as.numeric(Param.struc.1000[,3,14])-1)^2)), x=mean(correlation.1000[,10]), y=sqrt(mean((as.numeric(Param.struc.1000[,3,18])-1)^2)), lty=2, col=rgb(255,165,0,250,maxColorValue=255), lwd=2)

# green
segments(x0=mean(correlation.1000[,6]), y0=sqrt(mean((as.numeric(param.selec.1000[,3,2])-1)^2)), x=mean(correlation.1000[,7]), y=sqrt(mean((as.numeric(param.selec.1000[,3,6])-1)^2)), lty=2, col=rgb(0,139,139,200,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,7]), y0=sqrt(mean((as.numeric(param.selec.1000[,3,6])-1)^2)), x=mean(correlation.1000[,8]), y=sqrt(mean((as.numeric(param.selec.1000[,3,10])-1)^2)), lty=2, col=rgb(0,139,139,200,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,8]), y0=sqrt(mean((as.numeric(param.selec.1000[,3,10])-1)^2)), x=mean(correlation.1000[,9]), y=sqrt(mean((as.numeric(param.selec.1000[,3,14])-1)^2)), lty=2, col=rgb(0,139,139,200,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,9]), y0=sqrt(mean((as.numeric(param.selec.1000[,3,14])-1)^2)), x=mean(correlation.1000[,10]), y=sqrt(mean((as.numeric(param.selec.1000[,3,18])-1)^2)), lty=2, col=rgb(0,139,139,200,maxColorValue=255), lwd=2)

# purple
segments(x0=mean(correlation.1000[,6]), y0=sqrt(mean((as.numeric(param.selStruc.1000[,3,2])-1)^2)), x=mean(correlation.1000[,7]), y=sqrt(mean((as.numeric(param.selStruc.1000[,3,6])-1)^2)), lty=2, col=rgb(160,32,240,200,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,7]), y0=sqrt(mean((as.numeric(param.selStruc.1000[,3,6])-1)^2)), x=mean(correlation.1000[,8]), y=sqrt(mean((as.numeric(param.selStruc.1000[,3,10])-1)^2)), lty=2, col=rgb(160,32,240,200,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,8]), y0=sqrt(mean((as.numeric(param.selStruc.1000[,3,10])-1)^2)), x=mean(correlation.1000[,9]), y=sqrt(mean((as.numeric(param.selStruc.1000[,3,14])-1)^2)), lty=2, col=rgb(160,32,240,200,maxColorValue=255), lwd=2)

segments(x0=mean(correlation.1000[,9]), y0=sqrt(mean((as.numeric(param.selStruc.1000[,3,14])-1)^2)), x=mean(correlation.1000[,10]), y=sqrt(mean((as.numeric(param.selStruc.1000[,3,18])-1)^2)), lty=2, col=rgb(160,32,240,200,maxColorValue=255), lwd=2)



legend(.4,.4, legend=c("Stc", "HSl", "StSl"), col=c(rgb(255,165,0,250,maxColorValue=255),rgb(0,139,139,200,maxColorValue=255),rgb(160,32,240,200,maxColorValue=255)), cex=1.2, pch=16, bty="n")
